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2 . We develop a calculational scheme in Coulomb and temporal gauge that respects gauge 

^ ' invariance and is most easily applied to the infrared asymptotic region of QCD. It resembles 

I , the Dyson- Schwinger equations of Euclidean quantum field theory in Landau gauge, but 
Qh. 

^ • is 3- dimensional. A simple calculation yields a color-Coulomb potential that behaves at 



large R approximately like K;oui(-R) ~ o.2{d i)] spatial dimension 1 < d < 3. This 



^ ! is a linearly rising potential plus a rather weak dependence on d. 



1. Introduction 

There is a simple confinement scenario in Coulomb gauge [1], [2] which, in short, 
attributes confinement to the long range of the color-Coulomb potential, V"coui(-R)- This 
quantity is the instantaneous part of the 00-component of the dressed gluon propagator in 
minimal Coulomb gauge, ^ 

5fo-Doo(^, xo) = (fi'Ao(f , xo)gAl{0, 0)) = Vcovl\{\x\)S{xo) + (non - instantaneous), (1.1) 

and is given by [3] 



Here M{A) = —diDi{A) is the Faddeev- Popov operator, and the gauge-covariant derivative 
is defined by [Di{A)uj]'' = ^jCj" + got^^A^u". 

We present a calculation of Vcou\{R)- This quantity is of interest because: (i) It couples 
universally to color charge, (ii) Confinement of color-charge may be explained by the long 
range of this potential, (iii) It is a renormalization-group invariant and is independent of 
the cut-off and the renormalization mass [2]. (iv) A necessary condition for the Wilson 
potential V{R) to be confining is that Fcoui(-R) be confining [4], and if both potentials rise 
linearly at large i?, V{R) ~ aR and Fcoui(-R) ~ c^coui-R, then ctcouI > cr- (v) We wish to 
compare with a recent numerical determination [5] of V^oui(-R)) that does show a linear rise 
at large i?, with Ucoui ~ 3cr. 

Calculations in the Coulomb gauge have been pursued vigorously. For recent work and 
further references, see [6]. The present approach is distinguished by particular attention 
to gauge invariance, and its easiest application is to the infrared asymptotic limit of QCD. 

2. Temporal gauge and Coulomb gauge 

For simplicity we consider pure gluodynamics. In the temporal or Weyl gauge, Aq = 0, 
the wave functionals ^{A) depend on Af{x) for z = 1,2, 3. The color-electric field operator 
is represented by E^{x) — i5/6Al{x), and the hamiltonian by 



In this equation x represents a 3-vector, but everywhere else in this article 3-vectors are 
represented by x. 



V^ouiix - y)5'^' = { 9l[M-\A){-d')M-\A)]-l ). 



(1.2) 




(2.1) 
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where Bf = eijkidjA"^ + ^gQf°''"^A'jA'^), and the f"-^"^ are the structure constants of the 
Lie algebra of the SU(N) group. Wave functionals in temporal gauge are required to 
be gauge invariant '^{^A) = '^{A), where g{x) G SU{N) is a 3-dimensional local gauge 
transformation, and ^Ai = gQ^g~^dig + g~^Aig. These continuum equations have precise 
analogs in lattice gauge theory, where the Kogut-Suskind hamiltonian replaces the Weyl 
hamiltonian. 

Poincare invariance of the continuum theory is preserved because the hamiltonian 
density = + -B^) satisfies the Dirac-Schwinger equal-time commutation relation 

[T'\x),T''{y)] = -t[T'\x)+T'\y)]d,d{x-y) + S.T., (2.2) 

where T°* = ^eijk{EfB^ + B^E^) is the Poynting vector [7], and S.T. is the Schwinger 
term. 

Inner products in temporal gauge, (\E'i,\E'2) — N J dA ^1{A)^2{A), are divergent 
because of the gauge invariance of the wave-functionals. They may be made finite by 
using the Faddeev-Popov identity 

(^r^,^2)= f dA^' detM{A'') <ifl{A'')-^2{A''). (2.3) 

The integral extends over 3-dimensionally transverse configurations in the fundamental 
modular region A, which is a region free of Gribov copies. To be definite we suppose that 
we are in the minimal Coulomb gauge, which is obtained by minimizing F^ig) = 
with respect to gauge transformations g{x), so ||A|| < for all g{x) and all A in A. 

Wave-functionals in minimal Coulomb gauge \1/(A*'^) are the restriction of gauge- 
invariant wave functionals in temporal gauge "^{A) to the fundamental modular region A. 
Conversely every wave-functional in minimal Coulomb gauge has a unique gauge-invariant 
extension to temporal gauge. Every point A in the interior of A, is a unique absolute min- 
imum (modulo global gauge transformations), so the strict inequality holds ||A|| < 
(for all g{x) that is not a global gauge transformation). But every point Ai on the bound- 
ary dA of A is related by a local gauge transformation g{x) to some other point A2 = ^Ai 
also on 9 A, with which it is degenerate, \\Ai\\ = 11^211- This gauge transformation may 
be infinitesimal, A2 = Ai + eD{Ai)uj, where D{Ai)(ju is tangent to dA. Gauge-invar iance 
requires that the wave-functional in Coulomb gauge be identified at corresponding bound- 
ary points, ^^(^2) = ^(^1)5 or, for the infinitesimal case, that the wave functional satisfies 
(D(Ai)oj, f^Ui) = 0. This provides the boundary condition that is needed to make the 
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hamiltonian in Coulomb gauge well-defined and symmetric. The identification of boundary- 
points is often ignored because one does not know explicitly what the boundary of A is. 
But in general it would be a violation of gauge invariance to ignore this identification and 
take arbitrary wave-functionals in so-called physical coordinates which are the transverse 
configurations in A. In order not to make this error, in the present article we shall use 
wave functionals that are manifestly gauge-invariant. 

As an example we exhibit an approximate vacuum wave functional that is gauge 
invariant. The variation of the color-magnetic field is given by 

5Bt = e,,fei?f ^ [b5A)1, (2.4) 

which defines the hermitian operator D{A) that is the gauge-covariant curl. Consider the 
wave functional 

^ = exp(^-^ J d^x [(L>2)-i/25]a^. (2.5) 

The operator D{A) has small eigenvalues when acting on longitudinal fields, but the 
Bianchi identity DiBi = insures that the wave-functional is regular. We have jj^ 
-[D{D^)-y^B]i^, and 

where f{x) = i[(l)^)^/^]""(x, 2/)|jy=a;, and ~ means that derivatives with respect to 
(L*^)"^/^ are neglected. The first term will cancel the magnetic energy density, which 
is the most singular term, being the product of quantum fields at the same point. We have 
f{x) = e + u{x, A), where e = ^[{d'^y^'^]f-^{x,y)\y=x is a divergent constant, and u{x,A) 
is a gauge-invariant non-local functional that vanishes with A. The Schrodinger equation 
reads 

H^^ J d^x [e + u{x,A)]^, (2.7) 
and is violated by non-local terms only. 



3. Calculational scheme 

The vacuum wave functional ^'o(^) is positive, and we write ^'o(^) = exp[— S'(^)/2], 
where S{A) is manifestly gauge invariant. We assume that S{A) is either an approximate 
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expression, such as the one given above, or a trial expression. With gauge-invariant wave- 
functionals it is difficult to evaluate matrix elements by direct integration, and we shall 
borrow techniques from Euclidean quantum field theory. We have |\&o(^)P = exp[— S'(A)], 
and we define the generating functional of equal-time correlators. 



normalized to Z{0) = 1. This is precisely the formula for the partition function or generat- 
ing functional of 3-dimensional Euclidean gauge theory in the minimal Landau gauge, with 
the Yang- Mills Euclidean action 5'ym(^) replaced by some gauge-invariant action S{A). 
Only the transverse part of the source Jf (a;) contributes, and we take Jj to be identically 
transverse, 5j Jj = 0, and Jj = J*"^. 

We don't have an explicit expression for A, and we rely on the argument of [8] that fun- 
damental modular region A and the Gribov region Q have the same moments or correlators 
so we may integrate over instead of A, 



Whereas A is the set of absolute minima of the minimizing functional, the Gribov re- 
gion Q is the set of relative minima. The matrix of second derivatives of the minimizing 
functional is the Faddeev-Popov operator M{A). It is a non-negative matrix at a relative 
minimum, so the Gribov region O is the set of transverse configurations A^^ for which all 
eigenvalues Xn{A^^) of M(yl*'^) are non-negative, Xn{A^^) > 0. The interior of O consists 
of points A*"^ where all eigenvalues are strictly positive, Xn{A^^) > (apart from a trivial 
null eigenvector corresponding to global gauge transformations) . Its boundary d^l consists 
of points where M{A^^) has a non-trivial null-eigenvector, M{A^'^)(jj = 0, so Ai(^*^'') = 0, 
and all other eigenvalues are non-negative, Xn{A^^) > 0, for A^^ e dfl. 

We don't have an explicit expression for Q either, but we may exploit the fact 
that the integrand of (3.2) vanishes on dO, to derive the Dyson- Schwingcr (DS) equa- 
tions nevertheless [9]. Indeed the Faddeev-Popov determinent vanishes for A^^ G dQ, 




(3.1) 




(3.2) 



det M{A^') = Yln^n{A^') = 0. Thus the identity 




dA^' 



S 




(3.3) 



= 



5Af{x) 
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(— ) Z{J) = Jf{x)Z{J), (3.4) 



holds without a contribution from the boundary dO,. The set of DS equations in functional 
form, 

follow from this identity. Here T,(A^'') = S'(^*'') - trlnM(^*'') is the effective action. 
Because the integrand vanishes on dCt, the DS equations have the same form as they 
would if the integral (3.3) were extended to infinity. It is not necessary to know the 
boundary dQ explicitly, and the cut-off at dQ is implemented by imposing on the solution 
of the DS equations the natural positivity conditions that must be satisfied by the (equal- 
time) correlator {Ai{x)Aj{y)), by the ghost propagator G{x — y) = {goM~^{A^^){x,y)), 
and by the higher-order correlators. 

As in Euclidean quantum field theory, we rewrite (3.4) as a functional DS equation 
for W{J) = InZ(J), which is the analog of the free energy 

By Legendre transformation we convert this to a functional DS equation for the analog of 
the quantum effective action r{A^^') = {A^^ , J) — W{J), where Af''^{x) = 

5S / ... _ 5 \ . 5T 



SAf'^ix) y ^ SA^) ~ sAf^y ^^'^^ 

where V is the gluon propagator in the presence of the source, V~^{A^'^) — j^j^- The 
problem of evaluating correlators by direct functional integration has been replaced by 
the problem of solving the DS equations and, given the action S{A), we can, at least in 
principle, calculate all correlators by solving the DS equations for r(^*^'^). 

Suppose we take a trial expression S{A, ^) for the gauge- invariant action that depends 
on some unknown parameters ^. These parameters are determined by minimizing E{^) = 
(H) = (*o, iy*o)- To calculate E{^), write H = + Hm- We have Hm = lJ dPx B'^{x), 
where Bf{x) = Bf{x; A^^), and the magnetic energy is given by 

Err^iO = (Hr.) = 11^^"^ (^^ Jj) Z{J)\j^O. (3.7) 

To calculate the electric energy 

E,{0 = {He) = \j detM(A*'^) \E{x)^^\\ (3.8) 



we evaluate in temporal gauge, and then restrict to With ^'o(^) = 6xp[— 
we have Ef{x)^o = ^5if^ = -i^ti^', ^)*o, where S^{x; A) = ^^g^, which gives 



Ee(e) = ^/^^^ / dA'' detM(yl*'^) £^{x;A'')eM-SiA'')] 



(3.9) 



Now E{^) = -Ee(0 + -^m(0 has, in principle, been expressed in terms of the ^ parameters 
that appear in S{A,^). 



4. Infrared Ansatz 

This program may be difficult to carry out, especially if the action S{A) is non-local. 
However if our experience with DS equations with action SYuiA) is a reliable guide, then 
a remarkable simplification occurs in the infrared limit, as we now explain. 

The DS equations with Yang-Mills action were first solved, with due attention to 
the ghost contribution, in [10]. The subject is reviewed in [11], and recent results are 
reported in [12]. We will follow the method of [13], [9]. It was found in these investigations 
that the ghost contribution is the dominant one in the infrared. For example, in the DS 
equation for the gluon propagator, the leading contribution in the infrared is provided by 
the gluon loop. It was subsequently realized [8] that in the DS equation (3.6), with effective 
action E = ^ym — trlnM, one obtains the correct infrared asymptotic limit by setting 
S'ym = 0. Thus the infrared asymptotic limit is entirely determined by the Faddeev-Popov 
determinant detM(A*'^) and the cut-off at the Gribov horizon dQ. One might think that 
the functional integral with 5'ym = would diverge. However the DS equations are merely 
a technique for evaluating the functional integral, and since they give a finite result with 
'S'ym = 0, it appears that cut-off at the Gribov horizon makes the functional integral 
converge. 

Infrared Ansatz: We shall assume that in the present case also, the correct infrared 
limit is obtained by setting S{A) — 0. Moreover once one sets S{A) = 0, the present 
calculation reduces to the calculation of the infrared limit in Landau gauge, where one 
has 5'ym(^) = 0, and we may use directly the solution of [9], where d now represents the 
dimension of space instead of the dimension of space-time. 

We briefly outline how the solution was obtained in [9]. The crucial point is that a 
solution was sought for which the ghost propagator, G{k), is more singular than 1/A;^ at A; = 
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0. This property has been called the "horizon condition", and it triggers the confinement 
scenario in Coulomb gauge. The horizon condition holds because the Gribov region f2 is 
bounded in every direction and, in a space of high dimension such as configuration space, 
entropy favors population density near the Gribov horizon dfl. At the horizon M~^(A^^) 
is singular, and this enhances the ghost propagator G{x — y) — {go{M~^)xy)) at large 
separation or small k. The coupled DS equations for the gluon and ghost propators were 
solved in [13], and [9], taking the tree-level expression for the ghost-gluon vertex, and 
imposing the transversality condition diAi = on-shell. 



5. Calculation of color- Coulomb potential 

To calculate Fcoui(-R) we write (1.2) as 

Koui(a: - y)S-' = Jdh{ g^%x,z;A'^) {-d^^) g'\z,y;A'^) ), (5.1) 

and g^'^x^y^A^') = go[M-'^(A^')]fy. Its expectation-value, G{x-y)5''^ = {g''\x,y; A^')), 
is the ghost propagator. We separate the expectation-value of the product in (5.1) into 
disconnected and connected parts, 

Vconiix -y)- j d^'z G{x - z) {-dl) G{z - y) + V,,^{x - y), (5.2) 

which reads, upon fourier transformation, 

KouKA;) = k^G'i}^) + V;on(A;). (5.3) 

The infrared asymptotic form of the gluon and ghost propagators depends on two 
infrared critical exponents. 



^ ^ (5.4) 



G'^^k') = 



By equating like powers of momentum in either the gluon or the ghost DS equation, one 
obtains in either case the same relation an + 2aG = — (4 — d)/2. We use this equality to 
eliminate ao, in favor of a = ao- In the infrared limit, only the ghost loop contributes to 
the DS equation for the gluon which reads gives (^d^g)"^ — ^oict, d), where by eq. (A6) 
of [9] 

T (r. ^ r(2« + l-d/2) V\-a + d/2) 
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The only process that contributes to the DS equation for the ghost propagator is emission 
and absorption of a gluon. In the infrared hmit this gives (^d^^)""^ = /g(q!, d), where by 
eq. (A17) of [9] 

r Njd-l) TT r(2« + l) T{-a + d/2) 

' ' 2 (47r)rf/2 sin(7ra) r2(a + l) r(-2a + (i/2) r(a + l + (i/2)' ^ ' 

We must solve 

d) = Ig{o(, d) (5.7) 
to find the infrared critical exponent a — a{d). 



6. Discussion of solution 

We are interested in spatial dimension 1 < d < 3. The integral for Id{c(, d) converges 
for a in the interval j{d — 2) < a < and Id{ol^ d) is positive in this interval and diverges 
at the end-points. The integral for Icioi, d) converges for a in the interval < a < 1, and 
it diverges at the end-points. However iQ^a^d) changes sign in this interval at o; = 
and is positive only for < a < jd. Thus we look for solutions for a in the range 
max[0, l{d-2)]<a< \d. 

First take d in the interval 1 < d < 2. We have \{d — 2) < so we restrict our 
consideration to the interval < a < ^d. From the values at the end-points it follows 
that there are an odd number of solutions, and from numerical plots one sees that there 
is precisely one solution a{d) for 1 < d < 2. At d = 1, lG{a,d) vanishes because of the 
coefficient {d — 1) in (5.6). This coefficient occurs because the ghost emits and aborbs a 
gluon whose propagator contains a d-dimensional transverse projector, which vanishes at 
d = 1. (The coefficient d — 1 is an exact property of the DS equation of the ghost, and 
holds also in a more refined evaluation.) In contrast, lD{oc,d) is finite at cZ = 1 because 
a coefficient d — 1 has been factored out of both sides of the DS equation for the gluon 
propagator. So a{d) vanishes linearly with [d — 1), and from (5.7) one obtains in the 
limit (Z — > 1, 

a((Z) ^ ^((Z- 1) 0.20264 (cZ- 1). (6.1) 

TT 

In fact this formula fits a solution a{d) of (5.7) to about 2 % accuracy in the entire interval 
1 < d<A. At d = 2 the exact solution is given by 

«(2) = \, (6.2) 
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which differs from (6.1) by 1.3 %. The fitting formula, 

anid) = lid-l), (6.3) 

represents a solution ai{d) to (5.7) with about 1 % accuracy in the interval 1 < d < A. 

Now consider spatial dimension 2 < d < 3. We have < j{d — 2), so we seek a solution 
in the interval ^(c^ — 2) < o; < jd. From the values at the end-points, there are now an 
even number of solutions, and from numerical plots it appears that there are precisely two 
distinct real roots ai{d) and a2{d) for 2 < d < 3, except possibly at d = ~ 2.662 where 
there appears to be one root at q;(2.662) 0.33095, which is thus a crossing point of the 
two roots. At the crossing point (6.3) gives ct/i (2.662) = 0.33240, which is accurate to 
0.5 %. At d = 2 the two roots are given by ai{2) — i and 0^2(2) = 0. Only the root 
ai{2) = i matches the one root in the interval 1 < < 2, so it is the physical one, and 
the root ai(2) = is spurious. Thus for 2 < d < dc ^ 2.662, the larger root is the physical 
one. For values d > dc above the crossing point we do not know which of the two roots is 
physical. At d — 3 the equality (5.7) simplifies to ^'^^''(^^2a)(jX2a)' "''^ ~ ^' ^^^^ roots 

ai(3) ^ 0.3976; ^2(3) = ^. (6.4) 

The fitting formula (6.3) gives a/i(3) = 0.4, which agrees with the first root to about 1 %. 
At d = 4, there are two roots, ai(4) = ^ 0.5953, and a2(4) = 1. The fitting 

formula (6.3) gives q;/i(4) = 0.6, still accurate to 1 %. The fitting formula for the second 
root 

aMd) = ^{d-2) (6.5) 

is exact at d = 2, 3, and 4. The two fitting formulas cross at d = | fa 2.666 and 
a = i ^ 0.333. 

From the relation + 2aQ = ^(d — 4), and the fitting formula (6.3), we obtain the 
critical exponent of the gluon propagator, ar){d) ~ —^ + j^{d—l). This has a rather weak 
dependence on the spatial dimension d and gives a gluon propagator D{k) that vanishes 
at /c = for 1 < (i < 3. Thus the would-be transverse physical gluon does not appear in 
the spectrum. 

The critical exponent of the color-Coulomb potential is defined by 

CliW = (P)^' (6-6) 
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Suppose for simplicity that we neglect the connected term in (5.3) in the infrared asymp- 
totic limit, leaving for another occasion an evaluation of this term. Then we have in this 
limit 

Koul(^) = (^p^l+2aG ' ''^■'^^ 

and we obtain for the infrared critical exponant of the color-Coulomb potential ay = 2aG- 
The color-Coulomb potential is given at large R by Koui(-R) ~ ^^-d+2av _ jf ^ggg ^j^g 
simple fitting formula (6.3) for ao, one gets ay = ^{d — 1), and Fcoui(-R) ~ Ri^-^-'^('^-^)\ 
which is a linear potential plus a rather weak dependence on d for 1 < d < 3. For 
comparison we note that if instead aaid) = \ (d — 1), which is not so different from our 
solution, then one gets for the critical exponent of the gluon = — |, and of the color- 
Coulomb potential ay = ^{d — 1). This gives an exactly linear potential V^oui(-R) ~ -R, 
asymptotically at large R. 

The second solution at d = 3, namely 02(3) ~ |, yields V^coui(^) ~ l/^^i which 
gives an exactly a linearly rising color-Coulomb potential. However this success must be 
regarded as partly accidental because our solution does not give an exactly linear potential 
at d = 2, and a correct calculation should work for both d=2 and d = 3. 

The deviation from a linear potential should not be regarded as a failure of the ap- 
proach because our truncation scheme requires making an educated guess for the ghost- 
gluon vertex. We have chosen the tree-level vertex, but different choices give slightly 
different critical exponents [10]. So in our approach there is an inherent uncertainty in 
the critical exponent of the color-Coulomb potential. Moreover a correction to (6.7) may 
result from a more accurate evaluation of (5.3). However, granted these limitations, our 
results are at least qualitatively correct, and capture the essential features of a confining 
theory. We used the horizon condition which is the qualitative requirement that the ghost 
propagator be enhanced in the infrared compared to l/k^. The dynamics of the DS equa- 
tion then determine the infrared critical exponents of the ghost and gluon propagators 
and of the color-Coulomb potential. These are consistent with the confinement scenario in 
Coulomb gauge, which requires an infrared suppressed gluon propagator and a long-range 
color-Coulomb potential. They are also in at least qualitative agreement with numerical 
studies [14], [5]. The color-Coulomb potential we have obtained is confining and not far 
from linear for 1 < d < 3. 
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